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FINAL  TECHNICAL  REPORT  (Dated  8/28/01) 

AFOSR  GRANT  VF4920-97- 1-0497,  P00002 

(150-6756) 

“SHOCK  PROPAGATION  AND  SUPERSONIC  DRAG  IN 
LOW  TEMPERATURE  PLASMAS” 

SUPPORT  FOR: 

AFOSR-SPONSORED  WORKSHOP  ENTITLED: 
“Understanding  and  Control  of  Ionized  High-Speed  Flows” 
Princeton  University 
Princeton,  NJ  08544 
February  26-27, 1998 


ABSTRACT 

This  grant  supported  the  costs  of  an  Air  Force-sponsored  workshop,  “Understanding  and  Control 
of  Ionized  High-Speed  Flows,”  which  was  conducted  at  Princeton  University,  February  26-27, 
1998.  In  recent  years,  there  has  been  a  great  deal  of  interest  in  the  formation  mechanisms  and 
properties  of  air  plasmas.  This  interest  is,  in  part,  motivated  by  experiments  conducted  in  Russia 
and  in  the  United  States  which  indicate  shock  propagation  in  weakly  ionized  air  plasmas  is  at  a 
higher  velocity  than  would  be  predicted  by  presently  understood  models.  If  this  is,  indeed,  the 
case,  such  plasmas  could  be  used  for  supersonic/hypersonic  drag  reduction.  In  addition, 
atmospheric  plasmas  could  influence  flow  control  devices,  electromagnetic  attenuation,  and 
hypersonic  propulsion  systems.  As  a  consequence,  the  formation  of  such  plasmas  in  atmospheric 
pressure  environments,  and  the  study  of  the  properties  of  these  plasmas,  are  of  significant 
national  interest. 

DISCUSSION 

The  purpose  of  the  workshop  was  to  highlight  the  basic  science  issues  involved  with  controlling 
high-speed  flows  using  weakly-ionized  plasmas.  The  specific  objectives  of  the  two-day  meeting 
were  to  review  recent  experimental/theoretical  research  in  this  area,  and  to  exchange  ideas/ 
understanding  on  the  basic  flow  physics  of  the  problem,  and  to  discuss  future  research  directions. 
Princeton  University  was  a  natural  location  for  this  workshop  since  work  is  underway  here  to 
examine  both  drag  characteristics  and  electron  discharge  physics  associated  with  weakly  ionized 
plasmas. 

Sixty  individuals,  each  of  whom  is  considered  a  leading  authority  on  this  issue,  attended  the  two- 
day  workshop  (a  full  listing  of  the  participants  is  located  in  Appendix  A).  The  first  day  of  the 
workshop  consisted  of  four  sessions  on  an  agenda  of  twelve  speakers  and  covered  the  following 
topics: 
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Session  I:  “Introduction  to  the  Problem”  (presentations  by  B.  Ganguly  of  AFRL  and  A. 
Auslender  of  NASA  LaRc). 

Session  II:  “Aerodynamics  of  Ionized  Flow”  (presentations  by  R.  Miles  of  Princeton 
University,  J.  Shang  of  AFRL,  N.  Malmuth  and  V.  Soloviev  of  Rockwell  Science,  and  V. 
Subramanian  of  The  Ohio  State  University). 

Session  HI:  “Theory  and  Computation  of  Ionized  Flows”  (presentations  by  G.  Kamiadakis 
of  Brown  University,  K.  Powell  of  the  University  of  Michigan,  and  S.  Nazarenko  of  the 
University  of  Arizona). 

Session  IV:  “Plasma  Generation  and  Maintenance”  (Presentations  by  S.  Kuo  of  Polytechnic 
University,  J.  Kline  of  Research  Support  Instruments,  and  W.  Rich  of  The  Ohio  State 
University). 

The  second  day  of  the  workshop,  the  participants  were  divided  into  four  separate  working  groups 
for  further  discussion,  followed  by  debriefings  and  laboratory  tours  at  Princeton. 

Appendix  A  is  a  complete  record  of  the  two-day  workshop  and  includes  the  final  agenda  for  the 
meeting,  the  complete  list  of  participants  and  their  contact  information,  as  well  as  copies  of  each 
of  the  presentations  for  all  of  the  four  topic  areas  listed  above. 


Princeton  University 


School  of  Engineering  and  Applied  Science 
Department  of  Mechanical  ana  Aerospace  Engineering 

P.O.  Box  CN5263 

Princeton,  New  Jersey  08544-5263 
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AFOSR-SPONSORED  WORKSHOP: 


"UNDERSTANDING  AND  CONTROL  OF  IONIZED 
HIGH-SPEED  FLOWS" 


February  26-27, 1998 

MEETING  LOCATION:  Convocation  Room  (C217)  Engineering  Quadrangle 

Olden  Street,  Princeton,  NJ  08544 


THURSDAY,  FEBRUARY  26, 1998: 


7:30  AM 
8:00  AM 


8:30  AM 

9:30  AM 
10:00  AM 
10:30  AM 


11:30  AM 
12:00  PM 
1:30  PM 

2:30  PM 
3:00  PM 
3:30  PM 

4:30  PM 
5:00  PM 
6:00  PM 


SIGN-IN/COFFEE 

OPENING  REMARKS:  G.L.  Brown-Chairman,  MAE  Dept.,  Princeton 

R.  Miles-Princeton 

S.  Walker-AFOSR 

SESSION  I:  INTRODUCTION  TO  THE  EEOBLEM 

B.  Ganguly— AFRL,  A.  Auslender— NASA  LaRC 

DISCUSSION 

BREAK 

SESSION  n:  AERODYNAMICS  QE1QN1ZED.  EWE 

R.  Miles-Princeton,  J.  Shang-AFRL, 

N.  Malmuth/V.  Soloviev-Rockwell  Science 
V.  Subramanian-OSU 

DISCUSSION 

LUNCH 

SESSION  III:  THEORY  ANL>  COMPUTATION  QF  IONIZED.  FLOWS. 

G.  Kamiadakis-Brown  U.,  K  Powell-UM,  S.  Nazarenko-UA 

DISCUSSION 

BREAK 

SESSION  IV:  PLASMA  GENEEAIIONAm  MAINTENANCE 

S.  Kuo— Polytechnic  Univ.,  J.  Kline-RSI,  W.  Rich— OSU 

DISCUSSION 

FIRST  DAY  WRAP-UP 

RECEPTION,  Prospect  House 


FRIDAY,  FEBRUARY  27, 1998: 


7:30  AM  SIGN-IN/COFFEE 

8:00  AM  OPENING  REMARKS 

8: 15  AM  WORKING  GROUPS 

10:30  AM  BREAK 

11:00  AM  GROUP  DEBRIEFINGS 

12:00  PM  LUNCH/LABORATORY  TOUR 


ENGINEERING  QUADRANGLE 


A— WING  Chemical  Engineering 
B— WING  Electrical  Engineering 
C— WING  SEAS  Faculty  Lounge 
D— WING  Mechanical  <5 c  Aerospace 
Engineering 

E-WING  Civil  Engineering  and 
Operations  Research 

G— WING_  Energy  Research  Lab 
H— WING  Von  Neumann  Hall 
J-WING  POEM/EE/MAE 
COS  Computer  Science 


VI  Elevators 


O  Entrances/Exits 
•  Handicapped  Access 


J-WING 


D-WING 


Room  #3  (Rm.  G-201) 

Eisenhart  Conference  Room 


Room  #2  (Rm.  J-223) 

MAE  Faculty  Lounge 


Main  Meeting  (Room  C-217) 

Convocation  Room  ' — 


Room  #1  (Rm.  C-222) 

Dean's  Conference  Room 


BfflB'RA'RYi 


C-WING 


E-WING 


Main  Entrance 


Olden  Street 


Restrooms  are  located  in  the  comers 
of  the  building  where  2  wings  meet. 
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Examine  Electromagnetic  Manipulation  Of  Ionized 
Flowfields  To  Improve  Hypersonic  Vehicle 
Performance 


NEW  WORLD  VISTAS 
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Technical  exchange  on  flow  physics 
Discuss  future  directions 
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EXPERIMENTAL  OBSERVATIONS  OF  PLASMA-SHOCK 

INTERACTIONS _ 
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CURRENT  STATUS  QF  THE  RESEARCH  RESULTS 
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Shock  Wave  Propagation  in  a  Weakly  Ionized 

Low-pressure  Discharge 


A.  I.  Klimov  et  al.,  Sov.  Tech.  Phys.  Lett.  8,  192  (1982) 


Shock  Standoff  Comparison  to  Predictions 

(0.75  in.  diam,  30  torr) 


bQ  g  Li  t  i  i  «  i  #  I  »  «  ■  •  t  «  i  I  i  y  JJJLl±UXftjLLU  1X1  I  I  I  I  I  1  I  l-l- 

IJtU  -0.2  0.0  0.2  0.4  0.6 

Distance  ffom  Stagnation  Point,  inch 


MACH  NO.  DEPENDENT  SHOCK  WAVE  ATTENUATION  AND  DISPERSION 
VS.  GAS  PRESSURE  AND  CURRENT  DENSITY 
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KINETIC  TO  ELECTROSTATIC  ENERGY  TRANSFER 
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Schematic  of  Experiment 


Local  Velocity  in  Discharge  as  Function  of 
..  Discharge  Current 

Puls#  Energy  36Joule  Stored  Pulse  Energy  144Joule 

5  Torr  Nitrogen  i  20  Torr  Nitrogen  ^ 
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Discharge  Current  (mA)  Discharge  Current  (mA) 


SHOCK  WAVI  INTERACTIONS  WITH  WEAKLY  IONIZED 
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OUR  PLANNED  RESEARCH  FOR  FY98-99 


SUMMARY 
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Long  range  cooperative  interactions  may  be  important 


A  Comparison  of  the  Shock  Thickness  With  a  0.8  mm 
Diameter  Wire  ,  30  Torr  Argon  ,  Puiser  Voltage  9kV 


Argon  30  Torr,  lOkV  2(iF,  42.2cm 


50  Torr,  42.2cm 


600  620  640  660  680  700  720  740  760 
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Electrical  Schematic 


Deflection  Signals  Downstream  of  Discharge 

Pulse  Energy  81  Joule 
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AFOSR  /  High-Speed  Workshop 


NASA  Langley  Research  Center 


AFOSR  /  High-Speed  Workshop 
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Absorption  Ability  of  the  Microwave  Streamer  Discharge : 


The  Experimental  Data  about  a  Shock  Wave  Attenuation 

in  a  Microwave  Discharge 
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The  Acoustic  Interferometry  Study  of  the 
Microwave  Discharge  Region 
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NASA  Langley  Research  Center 


Plasma  Drag  Reduction  Physics 
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Plasma  Drag  Reduction  Physics 


Development  of  related  mechanism(s)  and  analysis 
of  the  relevant  “fluid”  equation  set  (OSU) 
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Shock  dispersion  experiments  and  plasma  generation 

(future  work) 


Anomalous”  Effects 
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B.  Supersonic  flow  in  weakly 
ionized  gas  around  a  spherical 
model  (ballistic  range 
experiments) 


Shock  Front  Velocity  (m/s) 


Shock  Propagation  Velocity  in  Weakly  Ionized  Gas 
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Charge  Transfer 
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EM  Field  at  DEL 


Electric  Potential  Distribution 
on  the  Centerline 


Shock  Structure 
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Fast  electrons  escape  DEL. 

Most  electrons  trapped  in  relaxation  zone 
Electrons  heated  by  inverse  bremstrahlung. 
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B.  Ionization  wave  generated  by  heavy  particles  (fast  ions  and 
neutrals)? 

C.  Onset  of  turbulence  by  angular  momentum  transfer  in  energy 
pooling  reactions? 
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Precursor"  Microwave  Plasma  at  Mach  6 


Precursor"  Microwave  Plasma  at  Mach  6 


Microwave  Breakdown  in  Air,  O,  and  N 


p  (mm  Hg) 

Cw  breakdown  fields  for  air  and  nitrogen  at  a  frequency  of  994  Mc/sec. 
(Characteristic  diffusion  length  A  =  2.65  cm) 

A.  D.  MacDonald  et.  al.  Phys.  Rev.  13£L  1841  (1963) 


Microwave  Breakdown  in  Air,  O,  and  N 


(Characteristic  diffusion  length  A  =  0.22  cm) 
A.  D.  MacDonald  et.  at.  Phys.  Rev.  130. 1841  (1963) 
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Program  Overview 

Shock  Wave  Dynamics  in  Weakly  Ionized  Plasmas 

at 

Laboratory  for  Modern  Fluid  Physics,  Florida  A&M  University  (FAMU),  Tallahassee,  FL 


QJ 

g  S 

«  QJ 

too  Y\  C5 

KT\ 

■O  js’M 

-a 

«.S.h 

Sa  S 

QJ  rj  43 

»*H  Tj 

o»  u 
CD  O  £ 

*  8* 

mh  d 
O  QJ 

*-<  ^  *d 
^  (Q  3 
T3  QJ  ** 
QJ 

J5  5  « 
«  K  *3 

•  pm  w  Cy 

S  d  3 
too.g  *3 

.S*  2 

»  JSO 

S  g;* 

a<  S 

g< 

5  to  § 

i-si 

Jat 

O  •  QJ 

•  n  cd 

rs 


44 

tQ 

QJ 

9h  £i 

3  0> 

^  04 

«£3  ^ 

cal  as 

£ 

•ts*« 

•  ^4 

r- 

.5 

£  S 

AS 

(0  d 

CA 

y-y 

QJ  <Q 

u 

i>  .. 

a 

£  .35  d 

£  d  3 

i«»*  «*■* 


'B.Mj*  g 

toO  v  u 

S3  O  0) 
oj3  2* 
.3  42  «  3 

9-4  A  CD 

r£>  **■'  O  S. 
•*2  cd  2 

S3  w  2  O 

is  o ^ 

cr*  S  «3h  qi 

S  d  y  £ 
o  2 

S'3  «  3 

,.«*-»  »^4 

V4  C3 

o  o  -S  “ 

S  CL,  qj  3 

^2  2  6 
C  d  *•  to 

rt  to  ^  JS 

k  ft  d? 
0J  ^  CL 

1  >•-  E 
.2  5  13  § 
a  ^  o-c 

i-Q  v>  G  _Q 

8*3  *ps 
5  o  .2  *d 

43  cn  CT1 

.2  C  «  2 

I  r -a  § 

£  S'd.B 

M’S  °  VI 

o»  ^  «  o 

o  S.rt’d 

^  «  tfl  r- 

o>  2s  ^  rt  . 

rl  toO  QJ 

S  m  43  qj 
7,  qj  S 
CU  ^  S  d 
2.2  CL  qj  ■ 

d  g  S  d 

i  ,,h  8  *° 


&P  13 

C3  ^ 


u  v-J 

o  <3 
W  ^VT3 

«  2  o 

“  S  E 

«2  E- 

£  o,2 

rj  y-» 

'g  §  c 

°  4-  QJ 
CD  CLi  n 

g-a  S. 

st-§  o 

cu  rt  cd 

ISSS 
aP  2  u  -3 

n  g  55 

5=5  is  8 

g  §5gj*i 

8-SSrfc  * 


duction  of  operational  aircraft. 


d  •- 

£)  QJ 

is  43 

5  d 

u 

M 
2  » 
M-l 

g  2 

O  3 
•»h  *3 
-*->  *3 
U  d, 

d  d 
H3  *■« 
Of  V 
Vi  V-i 

tS  ^ 

3  53 

01  CL, 


•a  | 

O  d 
Jd  *E 
w43 

MH  C3 

°  *3 

Q)  rr1 
u 

d 

QJ 

TJ 


D 

01  O 


c/5  8 

CJ  UH 

3  to 

°  6 
2  « 

QJ  f« 

CL|  rv 

5- 

&  i 

o>  53 

d  ^ 

g 
3  > 

O  o 

d  s 

,rH  QJ 

d  £ 

o  d 

*«s  ^ 

QJ  CJ 

i— < 

a  ® 
u  45 
u  co 

(M 
QJ  w 

>•6  » 
5  «  S 
SS® 
,y  „  — 

KJ  CL 

o  g_7 
8 

'g  £  £ 
>» « 
£33,3 

a 

£  «  £ 
’urr.  r 

•  ^  QJ 

«  |  § 
bis 

2  Sj  ^ 


■  CD  *tS 

s  §  J 

0:3  d 
ns  *53 

^>5  & 

f-  U  . 

M>S 

*P  53  ^ 

^  3  ° 

'S  «J3 
O  3  ’is 

«-a  g> 

■§•§  I 

lh  ^ 

“  ^  CO 

44  cr1  qj 
y  QJ.D 
o  d  *j 
•d  o  3 
w  c  5 

|  s  s 

•as  S 
D-Sfl 

SS" 

<3  v  J5 

i—i  cj 

V  3  »3 
*3  W  3S 

§  ”  s 

2«  -a 


3  6° 

CU  «  CO 
K  CO 

qj 

ixi  g  JS 

«_,  O  **"' 

=  V<  M-l 

3^0 
°  v«  > 

s^t 
s  s« 

^  S3  w 
g  **  T3 
s  gsS  aj 
•3  S?g>  3 

1  g 

a>  oi  J2  ^ 

b'SSo 

2  N  c-S 

CJ  CO 


a> 

u 

C  co 

QJ  O) 

EJ  .ts 

•el 

2  3 

cr 

u-( 

o  tS 
o 

£  CU 

y  to 

S  ® 

2  5 

■d  jy- 

g  « 

.s  *s, 

«  8 
43  55 
~  2 
m-i  y 
o  ta 

■»! 

§1 

•  *H  jo 

5  *d 

to  4-» 

0>4  ^ 

X  g 

01  o 

co 

43  o» 

o  QJ 

rt  ia 

to  <u 

6  3 

*b*a 

a-* 

C  *2 

qj  p 

d  cr 

.2  qj 

Vi  3 

«o 

Du  3 
X  , 

W  O 

b^ 

«  d 
d  *« 


A 
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AERODYNAMIC  APPLICATIONS 
OF  WEAKLY  IONIZED  PLASMAS 

R.  Miles,  S.  Macheret,  S.H.  Lam,  R  Efthimion,  L  Martinelli 
Princeton  University 


Drag  Reduction 
Plasma  Ramparts 

Flow  Control 

By  Localized  Energy  Addition 
By  J  x  B  forces 

Control  and  Intiation  of  Chemical  Processes 

Power  Extraction 

MHD  in  an  engine  inlet  (AJAX) 

MHD  in  ionized  flow 

Hypersonic  Wind  Tunnel  Drivers 

By  Energy  Addition  through  microwaves 
By  MHD  acceleration 

Suppression  of  shock  strength  and  sonic  boom 


SUSTAINING  A  PLASMA  IN  HIGH  SPEED  FLOW 


Important  for  Drag  reduction,  Plasma  Ramparts,  Flow  Control, 
Shock  Suppression,  MHD,  and  Combustion  Initiation  and  Control 


Discharge  may  be  blown  out 


Thermal  (equilibrium)  plasma  will  need  to  be  reinitiated  - 
a  thermal  plasma  cannot  propagate  at  supersonic  speed 


Nonequilibrium  plasma  will  require  high  electric  field  and  high  power 


A  filamentary  plasma  may  serve  the  purpose  and  will  require  less 
power  (a  small  portion  of  the  volume  is  ionized) 


Filaments  may  be  guided  by  laser  or  electron  beams 


A  high  power  microwave  may  be  the  preferred  source. 


SUSTAINING  MICROWAVE  DRIVEN 
PLASMAS  IN  ATMOSPHERIC  PRESSURE  AIR 


Thermal  Regime 

Requires  Tmin  >  4000K 

Needs  relatively  low  microwave  power  (kilowatts/cm2) 
the  power  offsets  the  thermal  losses 

Long  lifetime  (cooling  time) 


Nonequilibrium  Regime 

Requires  high  power  microwaves  to  overcome  the  electron 
recombination  and  attachment  rates 

Diffuse  plasma  (low  pressure) 

Low  temperature  neutral  species 


Filamentary  discharge  (high  pressure) 

High  temperature  localized  streamers 


FILAMENTATION  OF  ATMOSPHERIC 
MICROWAVE  DRIVEN  PLASMA 
Grachev,  Esakov,  Mishin  and  Khodataev 

Moscow  RadioTechnical  Institute 
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drag  reduction 

Experimental  Observations 


Shocks  Accelerate  Upon  Entry  Into  Plasmas  From  Neutral  Gas 


Dramatic  Changes  In  Structure: 

Shock  Front  Broadens 

Amplitude  Decreases 

Shock  Typically  Splits  Into  Two  or  Three 

Effect  Persists  for  Quite  a  Long  Time  (up  to  1  - 10  msec) 


Trans'vereenn^arirf13?  V(a‘K>US  Dlscharges  (Longitudinal  And 

Using  (  U'’  A‘SOn'  Hydr°gen'  HeIiU™'  Ete'> 


Ballistic  Range  Experiments  Seem  to  Imply  a 
Effective  Speed  of  Sound,  and  Show  Decrease 


Change  in  the 
in  Supersonic  Drag. 


Thermal  diaphragm 
Density  Contours 
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Length 

Planar  Cut  @(  0.83,-0.24,  0.25)(  0.83,  0.74,  0.25)(  1.82,-0.24,  0.25) 


Length 

Planar  Cut  @(  1.01,-0.28,  0.25)(  1.01,  1.26,  0.25)(  2.55,-0.28,  0.25) 


EXPERIMENTAL  CONCERNS 


Temperature 

In  all  the  experiments,  the  temperature  of  the 
neutral  gas  has  not  been  directly  measured. 


Temperature  Profile 

A  transverse  temperature  profile  or  a  curved 
thermal  boundary  will  lead  to  complex  shock 
structure.  This  may  occur  as  the  shock  passes 
through  the  hot  cathode  or  anode  fall  regions 
of  the  discharge.  (Well  away  from  the  centerline 
of  the  experiment) 

Models  indicate  the  apparent  shock  splitting 
and  the  density  rise  echo  experimental 
observations 


'Nonuniform  Temperature  Field 

Shock  splitting  has  been  seen  in  turbulent 
mixtures  of  gases.  If  the  plasma  has  hot  and  cold 
regions,  a  similar  shock  break-up  may  occur. 


CONCLUSIONS  ON  DRAG  REDUCTION 


All  the  experimental  phenomena  we  examined  are 
consistent  with  thermal  gradient  effects 


Ballistic  range  data  needs  to  be  studied  in  more  detail 


Accurate  temperature  and  density  measurements  need 
to  be  made 


There  may  be  some  effects  of  thermal  gradients,  unsteady 
phenomena,  vorticity,  and  random  or  flickering  heat  addition. 


PRINCETON  PROGRAM 


Development  of  new  diagnostics  for  temperature  and 
density  imaging 


New  World  Vistas 

Sustaining  microwave  driven  plasmas  in  a  Mach  3  flow 
Plasma  effects  on  the  bow  shock  in  a  Mach  3  flow. 
Study  of  basic  mechanisms  of  acoustic  wave 
interactions  and  shock  dynamics  in  stationary  plasmas 


Plasma  Ramparts 

Experimental  studies  of  the  control  of  plasma 
filamentation 

Kinetic  model  development  for  pulsed  ionized 
air  plasmas 

Modeling  nonequilibrium  molecular  processes 
High  power  microwave  driven  discharge  dynamics 
Stabilizing  and  extending  microwave  driven  plasmas 
with  magnetic  fields 


Mariah  II  (Radiatively  Driven  Wind  Tunnel  augmented  by  MHD) 
Sustaining  conductivity  in  atmospheric  pressure  air 
with  electron  beams 

Acceleration  and  deceleration  of  supersonic  flows  in 
MHD  channels 

Energy  coupling  into  high  pressure  air 
MHD  acceleration  with  filamentary  discharges 


PLASMA  TEMPERATURE  MEASUREMENT 
BY  FILTERED  RAYLEIGH  SCATTERING 


Rayleigh  Scattering 


Total  Scattering  is  Proportional  to  Density 

Linewidth  is  proportional  to  Jf 
Signal  is  weak 

Background  scattering  is  at  nearly  the  same 
frequency 

Filtered  Rayleigh  Scattering 

Eliminates  Background  Scattering 
Permits  Measurement  of  Linewidth 


f  Laser 


A 

b 

Rayleigh  scattering  +  background 

c 

filter  profile 

d 

imaged  signal  (no  background) 

FILTERED  RAYLEIGH  SCATTERING 
IMPLEMENTATION 


Source 


Filter 


Camera 


Tunable  Ultraviolet,  Injection  Locked  Laser 
Ti: sapphire  with  “Ramp  and  Lock”  Technolo 
(Schwartz  Electro-Optics) 


Atomic  Mercury  Vapor 

2”  Diameter  Cell  Held  at  40°C  with  a 
hot  water  bath  (vapor  pressure  1  Torr) 


UV  Sensitive,  Intensified  CCD  (Princeton 
Instruments) 


SEQUENTIAL  RAYLEIGH  IMAGES  OF  MACH  2.5 
SHOCKWAVE  /  BOUNDARY  LAYER  INTERACTION 

AT  A  14°  WEDGE 


Framing  rate  is  500,000  images  per  second 


Mariah  II 


No  wind  tunnel  can  at  present  be  run  to  accurately 
simulate  flight  conditiona  at  higher  than  Mach  8 

240K,  (120,000  ft)  requires  2800K,  8  atm  in  plenum 
NO  contamination 

Low  static  T  or  low  static  P  (enthalpy 
too  low  or  entropy  too  high) 

Short  time  operation 

Radiatively  Driven  concept  adds  energy  downstream  of 
the  throat  in  the  supersonic  section 

Static  temperature  stays  low 
Low  temperature  but  high  pressure  in  the  plenum 
Concept  is  driven  by  lasers  or  electron  beams 
Demonstrated  with  10KW  laser  in  December  at 
Wright  Labs 

MHD  section  is  to  be  added  to  further  accelerate  the 
flow  after  the  Radiative  section 

Reduces  the  front  end  pressure  requirement 
Extends  the  envelop  of  the  tunnel 
Requires  electron  beam  sustained 
conductivity 

Must  operate  at  close  to  room  temperature 
and  at  pressures  on  the  order  of  .1  atm. 


Energy  Addition  by  Lasers 


RDHWT  Review  April  23  -  24,^1997 


Figure  1.  Diagram  of  10  kw  laser  energy  addition  experiment 


X  (m) 


0)  00  00 
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Static  1  inch  from  exit 
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Time  (s) 


Flow  Around  Probe 


Results  verified  by  Pitot  Measurements 


RDHWT  Review  April  23  -  24,  1997 
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Thermodynamic  limit  of  RDHWT: 

Ah  =  Ta& s 

Ta  can’t  be  extremely  high  (erosion,  chemistry); 
minimum  s  is  defined  by  incompressibility; 
maximum  s  is  defined  by  test  section  requirements. 

THERMODYNAMIC  ADVANTAGE  OF  MHD 

Total  enthalpy  augmentation: 

Ah.  =  — —  TAs 

1  K- 1 

where  K  =  E/(uB)  =  (Energy  added)/(Push  work) 
The  best  thermodynamically: 

when  A ht  can  be  large  even  for  small  T  As 

(non-entropy-generating  energy  addition). 

However,  when  JT— >1,  j  -  <j  (E-uB)  ->0, 
and  to  add 

Aht  =  cK(K-])uB2L/p, 

very  large  conductivity,  a,  and/or  length,  L,  would  be 
required 


AJAX 

Proposed  by  Russians  in  1990 
(Vladimir  Freistadt  - 

State  Hypersonic  Systems  Scientific  Research  Enterprise 

St.  Petersburg) 


MHD  converter  in  the  inlet  connected  to  MHD  accelerator  in  the 
nozzle  to  by-pass  free  stream  kinetic  energy  and  reduce  the  entropy 
rise  in  the  inlet. 

Reduces  the  flow  to  subsonic  in  the  combustor 
MHD  conditions  almost  the  same  as  Mariah  II 
Requires  some  method  to  sustain  conductivity 


Reduction  in  vehicle  drag  by  cold  plasma  injection  at  the  nose 


Increase  of  combustion  volume  and  efficiency  within  the  engine  by 
plasma  injection  or  injection  of  materials  ahead  of  the  fuel  injectors 


Reforming  of  kerosene  and  water  into  hydrogen  and  CO  for  fuel. 
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Noseboom  Plasma  Flow  Control 
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-  hypersonic  applications 
P4  study  useful  since 
experimental  dataset 
available  to  test  models 


Problems 


Problem  3a  -  (F-15) 
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1-D  shock  propagating  into  a  plasma 


Perfect  Gas  Shock-Layer  Problem 


6  =  viscous  shock  layer  tliickness 
Vf  -  Vt  s  shock  strength 
Y  s  specific  heat  ratio 


Scaling  of  Shock  Attenuation  with 
Thickness  in  a  Plasma 


Science  Center 


Drag  Reduction  Mechanisms 
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Flow  Configurations 


Aerospike  Effect  on  Pressures 


Peak  attenuates,  with  increasing  spike  length  reducing  drag 

Needs  to  be  traded  against  high  heat  transfer  near  peak 

Plasma  heating  can  affect  reattachment  and  therefore  pressure  peak 
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Comparison  of  Predicted  Drags  With  Russian 
Experiments 


Forebody  drag  of  a  cone-cylinder  with  counterflow  jet  in 
supersonic  flow _ _ _ _ 


Science  Center 


Forebody  drag  of  a  cone-cylinder  with  counterflow  jet  in 
supersonic  flow _ •  _ _ 
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vorticity,  multiple  secondary  shocks  and  instabilities 


Forebody  drag  of  a  cone-cylinder  with  counterflow  jet  in 
supersonic  flow _ _ _ _ _ 


Science  Center 


Figure  7:  Forebody  drag  coefficient  as  a  function  of  p4/pco  for  different  values  of  r4/T’00  as 
follows:  crosses:  1,  triangles:3,  diamonds:6,  asterisks:10.  The  forebody  drag  coefficients  of  the 
sharp  cone  and  the  cone  frustum  are  shown  as  lines  for  reference.  =  2,  7  =  1.4 


Figure  8:  The  minimum  value  of  the  forebody  drag  coefficient  as  a  function  of  Ta/T^. 
Mao  =  2,  7  =  1.4 
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Propagation  of  1-D  Shock  Waves  Through  Weakly  Ionized 

Plasmas 

V.R.Soloviev,  V.M.Krivtsov,  A.M.Konchakov 

Moscow  Institute  of  Physics  and  Technology,  Russia 

N.Malmuth 

Rockwell  Science  Center,  CA,  USA 

Main  topics  of  the  model 

•  Hugoniot  relations  on  a  shock  wave  (SW)  in  a  weakly  ionized  plasma 

•  Estimations  of  the  energy  stored  in  different  excited  states  for  glow  discharge 
plasma 

•  The  existence  of  high  situated  metastable  states  in  Ar  and  N2 ,  which  serve  as 
a  ‘new  ground  state’  for  higher  states  excitation  by  electron  impact 

•  Highly  excited  atomic  and/or  molecular  states  store  sufficient  energy  to 
produce  SW  acceleration  and  attenuation 

•  Inclusion  of  highly  excited  atom  energy  release  behind  the  SW  front 
dramatically  improves  agreement  of  theory  with  experiment  for  1-D  shock 
penetrating  into  a  plasma  and  hypersonic  sphere  entering  a  plasma 


Hugoniot  relations  on  a  shock  wave  (SW)  in  a  weakly  ionized 

plasma 


A  set  of  gas  mass,  momentum  and  energy  conservation  equations  in  SW 

coordinates 

£p+£^)=0> 

dt  d$ 


d(p(u  +  D))  d  ,  ,  .  d 

■  . ■+—(M"  +  -D)+p)  =  — 

dt  dd,  dd, 


4  dit 
—  u — 
3  d£ 


d 

[7a+o»+£)!] 

d 

i  _ 

Lf/,+(,,+Z))2T 

d 

'  A  dh  \ 

dt 

[T  2  J 

P 

^  \  2  j 

1 

& 

4  d 


2  r 


(u  +  D) 


+f'(u  +  D)  +  Or 


h 


=  Cja  = 

P  <* 


7  P 
Y-\py 


p,  p,  h  are  density,  pressure  and  enthalpy  of  the  gas;  v  is  its  velocity; 

A,  p  are  thermal  conductivity  and  dynamic  viscosity  coefficients; 

Cp,  y  are  specific  heat  at  a  constant  pressure  and  the  ratio  of  specific  heats; 
D  is  the  SW  velocity; 

The /*  and  Or  terms  describe  the  plasma  influence  on  the  gas. 


Hugoniot  relations 


pxu\  +  px  =  p0D2  (1  -eul)+ p0 , 

/h+i-=h,(U-st)+^-0-2s,2), 


where 


€ 


ul 


Su2  ~ 


1 

PoD 3 


J/^> 


£„  = 


1 

PoDh0 


\qm 


•  To  produce  a  plasma  effect  on  SW  propagation  it  is  necessary  £ui,  Su2,  or  q,  to 
be  of  the  order  of  1. 


The  typical  cold  glow  discharge  plasma  parameters  in 
SW  -  discharge  interaction  experiments 


gas  pressure  (gas  density) 
gas  temperature 
current  density 
electric  field  strength 
E/N  value 

discharge  specific  power 
electron  temperature 
electron  density 

metastable  atom  (or  molecule)  density 


p  «  3  -  30  Torr  ( N  «  1017  - 1018  cm*3 ) 
T  «  300  K 
j « 1  - 100  mA/cm2 
E  *  100  V/cm 
(2  -10)  xlO*16  Vxcm2 
W«0.1  - 10  W/cm3 
Te «  1  -  3  eV 
ne  «  1010  - 1012  cm*3 
Nm*10n-1013  cm*3 


The  energy  stored  in  different  excited  states 


•  freestream  gas  enthalpy 


ho=NcpT~1016  eV/cm3 


•  ionization  iijEi  ~  1012  eV/cm3,  8h=  n,Ei/  NcpT  ~  10"4 

•  metastable  atoms  (or  molecules)  NmEm~  1014  eV/cm3,  8h=  NmEm  /NcpT~10'2 

•  specific  discharge  energy  supply 

during  the  time  of  unit  length  SW  pass  £h=  Wx  lcm/D  NcpT~10*2  - 10*1 

•  the  vibrational  energy  is  shown  to  be  unsuitable  to  explain  the  SW  anomalious 
propagation  effects  because  of  large  vibrational-translational  relaxation  time 
tvt~10'3s 


The  scheme  of  Ar  and  No  electronic  excited  states 


A-+(15.76eV) 


N/15.5&V) 


4p(B2) 
4(11.8)  — 


04(11.0) 

HT|(735) 


a'rp.54) 


A^617) 


Ar(4s)  and  N2(A3E)  are  the  metastable  states 


The  Ar  metastable  state  quenching  and  building  up 


Process 

Ar(4s)  +  2Ar  -»  Ar*2  +  Ar 
Ar(4s)  +  e  -»  Ar  +  e  +  1 1.8ev 
Ar(4s)  +  N2  ->  Ar  +  N2(C,  B, ...) 
Ar(4s)  +Ar(4s)  -»  Ar+  +  Ar  +  e 
Ar(4p)  +  Ar  —>  Ar(4s)  +  Ar 
Ar(4p)  -»  Ar(4s)  +  hv 


Rate  constant  Quenching  time  for 

p=10  Torr,  ne=10ncm'3 
and  1%  N2  diluent  in  Ar 
10'32  cm6/s  0.8  ms 

10‘8  cm3/s  1ms 

3.6xl0'ncm3/s  0.01ms 

10'9  cm3/s  lmsx(1012/n(4s)[cm‘3]) 

~  10'11-  10‘12cm3/s  0.3  -  3mcs 

3.2xl07  s'1  0.03mcs 


The  N2(A3S)  metastable  state  quenching  and  building  up 


N2(A3S)  +  Ar  — >•  N2(X)  +  Ar 
N2(A3Z)  +  N2(X)  ->  2  N2(X) 

N2(A3X)  +  02  ->  N2(X)  +  02 

N2(A3X)  +  N2(X,v)  N2(B3n)  +  N2(X) 

N2(A3S)  +  N2(A3Z)  N2(B,C,...)  +  N2(X) 

N2(B3n)  +  N2(X)  N2(A3E)  +  N2(X) 

N2(B3ri) -»N2(A3S)  +  hv 
N2(C3n)  -»  N2(A,B)  +  hv 


<4x  1 0‘14cm3/s  >  0. 1ms 

<4.5xl0‘17cm3/s  >100ms 

3.6x1 0*12cm3/s  4mcs 

(for  air  composition) 

~  10“10cm3/s 

10'9cm3/s  lmsx(1012/n(A)[cm'3]) 

2xl0*ncm3/s  0.15mcs 

1.5x10s  s'1  6.5mcs 

2.7xl07  s*1  0.037mcs 


The  discharge  fractional  power  input  into  different 
excited  states  for  Ar,  N0  and  Ar-No  mixture 


J9JSUBJ}  J9M0d  |BUO|}OBJ  J 


E/N,  10-18V.cm 


The  mean  electron  energy  (electron  temperature)  and 
electron  drift  velocity  versus  E/N  value  for  Ar,  Ar+O.OIN 


ca  io  cm  ^  o 


A9  ‘A6J0U0  UOJP0|0  UB01AI 


The  kinetic  model  of  energy  storing  in  highly  excited  states 

i _ 


AE^ 
highly 

excited  states 
energy  store 

N* 


energy 
from  the 
discharge 


N** 


a  (to  heat 
tand  radiation) 

Nm  (quenching) 


The  system  of  balance  equations  for  excited  states 


^  ~  N.  +  — )  -  K«.Nm  +  N\k,N  +  Ak$) 

**  hm  Tq 

m' 


a 


=  KneNm-N\kqN  +  ARe) 


n /  &  E  —E  a  F 

—  exp( — 2 - £-),  for  n>n*  k^ie»k^[ — 

«=«*  8„-  Te  Tt 


^  2Ry 


A£\  .x,  «  n*«3 


n’"+1  «3  * 


\2Ry  KN 
T  kcN. 


*4-5,  ErE*  0.5-leV 


Ry  =13.6eV  is  the  Rydberg  constant 


W = jE  -  is  the  discharge  specific  power; 

tj~1  -  is  the  discharge  fractional  power  input  into  excitation; 

Em,  En  -  are  the  energies  of  the  different  states  excitation  from  the  ground  state; 
kp  -  is  the  metastable  atom  Penning  process  rate  constant; 
xq  -  is  the  metastable  state  quenching  time  in  other  processes; 
ke  -  is  the  electron  excitation  rate  constant  from  metastable  to  N*  states; 
kq  -  is  the  gas  quenching  rate  constant  for  N*  states; 

Ar-  is  the  radiative  decay  rate  for  N*  states; 

8  -  is  the  probability  of  radiation  trapping; 

AEj  -  is  the  decrease  of  ionization  potential  due  to  excited  atom  electron  -  ion  interaction; 


0  = 


1 


For  the  tube  radius  R  0  /  Ar  n 

D^7UG  mR 

For  Ar(4p)  -»  Ar(4s)  transitions  0~3xl  O'2  (1012/Nm)0  :)  followed  by 

AR0«kqN 


The  steady  state  values  of  Nm  and  N* : 


N 


m 


-—A  *2.5x10 B y[w\}Fum*], cm 

\kPEm 


-3 


kenc 


N  *  3.5x10 


12 


/w 


■Jw[W/cm3] 

/?[7b/r] 


,  cm 


-3 


The  equilibrium  highly  excited  hydrogen-like  states 
population 


-  the  upper  limit  of  principle  quantum  number  for 


hydrogen-like  states; 


nj&  100  forNi«10ncm'3 

....  „.2  n)  E,-E\ 

vexp( — ~ 1(,4N* 

g*  ~30  is  die  degeneracy  of  N*  state 


The  energy  stored  in  highly  excited  states 


s**  &  (E*  -  exp(-^y^) 
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The  energy  released  to  heat  behind  the  shock  wave 
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Inclusion  of  Plasma  Effects  Dramatically 
Improves  Agreement  of  Theory  with  Experiment 
for  1-D  Shock  Penetrating  Into  a  Plasma 
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Inclusion  of  Plasma  Effects  Dramatically 
mproves  Agreement  of  Theory  with  Experiment 
for  Hypersonic  Sphere  Entering  a  Plasma 
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Conclusions 


•  The  mechanism  of  shock  wave  (SW)  acceleration  and 
attenuation  due  to  interaction  with  a  weakly  ionized  plasma 
really  exists 

•  Preliminary  calculations  of  1-D  SW  speed  in  plasma  are  in  a 
good  agreement  with  experimental  results  when  the  highly 
excited  states  energy  release  behind  the  SW  front  is  taken 
into  account 

•  SW  acceleration  is  accompanied  by  decrease  of  density  and 
velocity  jump  with  a  small  increase  of  static  pressure  jump 
through  the  SW  front 

•  The  foregoing  results  in  a  decrease  of  dynamic  pressure 
jump  through  the  S W  front 

•  The  creation  of  a  discharge  plasma  IN  FRONT  OF  the  shock 
wave  courses  the  SW  attenuation  and  drag  reduction  for 
supersonic  and  hypersonic  vehicles 

•  More  detailed  kinetic  analysis  is  necessary  to  obtain  more 
precise  results 
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Source  terms  in  continuity  equation 
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shock  nearly  recovers  its  original  strength 
well  downstream  of  the  discharge. 


explain  some  experimental  observations, 
but  cannot  explain  others  (such  as  shock 
recovery  distance). 
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Flow  Control  via  MHD  and  EMHD 

George  Em  Karniadakis 
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High-Speed  Flow  Features 


Strong  shock  wave  effects/interactions 
-  thin  layers,  small  detachment  distance,  etropy  layer 
Strong  viscous  effects/interactions 


-  increased  drag,  heating  at  leading  edge 

Transitional  flow 

-  Rec  ~  108  vs.  105  (subsonic) 

-  Cj  and  Ch  are  3x  laminar 

-  Reec  ~  100 Me 

Nonequilibrium  ’*  '  MuUt  -  PaJ-<L  PhysicS 


-  thermal,  chemical 

-  introduction  of  length  scale 

-  multi-rate  physics 

Ionization 


-  4000  K  to  6000  K:  NO  -+  NO +  +  e"  (mild) 

-  >  9000  K :N->N+  +  e~ 

O.-*  0+  +  e- 


-  Above  20kM  +  10  km /s 

-  External  B,  E.  frequencies 

Multiple  boundary  layers 


-  hydrodynamics,  thermal,  species 

-  electric  conductivity,  etropy  layers 


Non-Equilibrium  Plasma  Boundary  Layer 

Effects  on  Te,  a 

atmospheric-pressure  potassium-seeded  arpon  plasma 


Normalized  Temperature 


Normalized  Conductivity 


Ref:  Brown  &  Mitchner  (1971) 

"...  Even  in  plasma  which  is  strongly  collision  dominated,  significant  nonequilibrium 
may  be  present  in  the  boundary  layer  adjacent  to  a  cooled  solid  surface ...” 


Note:  Joule  heating  due  to  E  could  further  increase  conductivity. 


Electrical  Conductivity 


•  In 


Note:  f 


ET 


Electric  Field  vs.  Ionization 


•onstant  E,  B  fields  tensor  conductivity 

o-||  =^(ne,mc,i/e/f,/°) 

Weakly-ionized  j  crx  =  T(ne,  mc,  veH,  /»,  ue) 

.  ~  J’iP'd  me,  Veffi  /*\  ^e) 

If5  =  0=J>-CrjL  =  (7||  =  (7  t 

If  E  ~  =»  dependence  of  <r  on  a;  }  f«t  ^  tM*-  ^  , 

9  =  » -E±>  T#,  me,  ) 


Lorentz  Force:  Numerical  Modeling 


•  Governing  Equations  for  B  and  E 
«  Boundary  Conditions  -  The  Motz  problem 
Example  1:  Electromagnetic  Hies 

Dirichlet:  FL  ~  Neumann:  F±  ~  F\\ 
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Example  2:  Alternate  stripes  of  electrodes/magnets  Fy 


~  e 


Dirichlet  Neumann  D-N  Mixed 


Lorentz  Force 


MHD  Assumptions  ( tc  >  upl) 

PcU 

-  Convection  current  neglected:  — —  <C  1 


-  Force  due  to  net  charge  neglected: 


PcE 


\J  x  B\ 


<  1 


Then: 


Fl  —  J  x  B\  Eq  =  external 
J  —  o{E  *f-  v.  x  B)  “h  o' Eq 


MHD-Control 

•  Apply  B 

•  E  =  — V$ 

•  V2  $  =  V-(wxB) 
oB2L 


•  ^induced  ~ 


pU 


EMHD-Control 


•4 

•  Apply  B  and  E 

•  J  ~  oEq 


•  I ; 


applied 


oEqBL 

pU2 


where  Rem  =  poUL  <  1  in  lab  scale  ^ac«|Uc*  i\4nca<£ 

•  But  Ohm’s  law  for  partially  ionized  gas: 

J  +  PeJ  xb  +  sb  x  (J  x  b)  =  a(J5  +  «xB)  +  ——'Vpe 
a.  ene 

J  “Generalized  Ohm’s  Law” 

Un-slip 


Discontinuous  Galerkin  Method  -  MekT ar 


dynamical  Direct  Numerical  Simulation  (dDNS) 

MHD,  Compressible  &  Incompressible  flows 

Full  3D  Configurations  -  Arbitrary  Geometric  Complexity 

Unstructured/Hybrids  Grids 


High-Order/Spectral  Accuracy 

Conservative  Formulation 

h-refinement  for  shocks/discontinuities 

ALE  Algorithm  for  Moving  Rotating  Subdomains 
METIS/MPI-based  Parallelization 


+  PE  *+  S 
2S. 
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W4CH2./S 
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,  Sp«M/Ap  •«  tfdurcd 
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H-P  Adaptive  Refinement:  Viscous  Supersonic  Flows 


Figure  1: 


Discontinuous  Galerkin  Simulation:  Density  contours 


Figure  2:  Unstructured  grid  and  variable-order  per  element 
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PhD  Thesis:  T.  Warburton 
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MHD/EMHD-Simulations 


I.  Laminar  &  Turbulent  Wakes 

£  -  vortex  street  suppression 

-  phase/frequency  control 

H.  Wall  Turbulence 

-  High  cr,  spanwise  forcing/MHD 

-  Low  cr,  steamwise  forcing 
f  -  Low  a,  “normal”  forcing 

Method:  Parallel  DNS  -  The  Afen'T ctr  Code 


Experiments:  Conflicting  results 


j*  FZR  8c  TU  Dresden 


T.  Weier,  O.  Posdziech 


EMHD  control  of  the  flow  around  a 

circular  cylinder 


L,  . .  - .  . . .  . .  . . . .  .  - - - - - - - - - - ; 

Dimensionless  numbers: 

Re=^d 

V 

-  Reynolds  number 

II 

•S’ 

ft. 

-  Interaction  Parameter 

P  voo 

Air  Force  Workshop 


Princeton  NJ 


Feb.  25,  1998 


Tt)R80t.e**T  AXAfic -  f jr 


Cylinder,  Re  =  500 


M-OA/ofa  /VI*  C.  J 


Active  Control:  Lorentz  Body  Force 


•*  Pk-i'  2  - *-  Electric  Field 

)  - ►  Magnetic  Field 

»  K  4  H  - ►  Lorentz  Force 


P  o£ 

4  TOS.R 


Henoch  &  Stace,  1995 


EMHD:  Shear  Stress  History  -  Effects  of  BC 
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NEAR- WALL  TURBULENCE  and  LAMB  VECTOR 


Instantaneous  distribution  of  turbulent  charge  (ch2). 


Figure  1:  Dipole  structure  of  the  divergence  of  the  Lamb  vector:  turbulent  charge 


z 


Figure  2:  Footprint  of  the  turbulent  charge:  LEFT  -  EMHD  on;  RIGHT:  EMHD  off 


Summary 


Explore  MHD,  EMHD  and  EHD 

Strong  electrical  conductivity  non-uniformity 

Selection  of  physical  model  is  application-dependent 

-  degree  of  ionization 

Transitional  Flow  -  No  ad  hoc  models 
Numerical  Model  (A/* ck'T otr  ): 

-  High  order 

-  Hybrid  grids 


Benchmark  Experiments 
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Right  Boundary:  p  =  o.i;  d  =  “V";  -a  =  v  =  o 
(top  and  bottom  B.C's  are  “reflecting  slip1* ) 

Initially,  th.es e  two  states 
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SIMULATION  METHOD; 

(of  the  2D  Euler  eqiations) 

PHYSICAL  ASPECT  RATIO  (x:y)  =  100:1 

#  of  gridpoints  in  x  direction  =  400 

#  of  gridpoints  in  y  direction  =  20 

#  of  ghostpoints  =3  (on  each  side) 


1 :  Given  the  fluid  parameters 
at  each  gridpoint,  find  the 
eigenvectors  to  Roe's  matrices. 

2 :  Use  a  5th.  order  weighted  ENO 
scheme  to  find  t  he  flux 
eigenvectors  at  half  gridpoints  • 


3:  Use  a  3rd  order  Runge  Kutta  routine 
to  propagate  the  half  gridpoint 
fluxes  and  find  the  new  fluid 
parameters  at  th  e  real  gridpoints  • 


4:  Repeat. 


PcoCt $U-*>  ^ro^lr  ot  "VUfl-  SUock 


/4JOO|9A  x 


X^pipoA 


Qjnss^Jd 


X}IOO  |3A  X 


sjnss&jd 


K}isus>p 


A}I00|3A  X 


ajnss&Jd 


i  ■  ■  ■  I  ■  ■  ■  ■  ■  ■  ■ _ l  ■  *  ■  n  0 


^  OM  O  CJ4 

/v}IQi:|JOA 


/qisusp 


002 


003 


001, 


0 


A}ic>o|&a-x  6ad 


.i L 


i L 


L 


tilt 


002 


003 


001. 


o  o  o 
A}!oo|sa— x  Bad 


00£ 


003 


001 


J  0 


I  ■  <  ■  I _ I _ J _ L _ I  I  0 


00£ 


003 


001 


djnssojd  Bad 


X^isusp  6 ad 


RSI  Further  Investigation  of  Large  Volume  PIA 

Research  (Persistent  Ionizationin  Air)  Plasmas  at 

SuPPort  Atmospheric  Pressure 
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Plasma  Torches  and  Their  Demonstration 


Spencer  P.  Kuo  and  Edward  Koretzky 


Department  of  Electrical  Engineering 
Polytechnic  University 


-*Work  supported  by  the  AFOSR 
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(3  x  3)  of  plasma  torches.  Tung-  plasma  torch, 
sten  electrodes  are  used.  Each  elec¬ 
trode  pair  is  facilitated  with  an  air 
jet  nozzle. 


oscilloscope 


03 


0.04 
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Plasma  Torch 
1  ATM 


Pressure: 

Density: 

Temperature: 

Dimensions: 

Volume: 

Power  Consumption: 

Peak  Voltage  and  Current: 


1.5  x  1013  (deduced  from  mea¬ 
surements  based  on  two  inde¬ 
pendent  methods) 

1200° K 

cylindrical  shape  having  1  cm 
radius  and  7  cm  length 

20  cm3/torch 

600  W/torch 

2.8  kV  and  2  A 


Microwave  Plasma  Layers 

Pressure:  20  torr 

Density:  1011  cm-3 

Microwave  Power:  1  MW  peak,  1  fis  pulse,  40  rps 


-  The  Ohio  State  University  jj - 

Plasma  Ramparts  Using 
Metastable  Molecules 

MURI  Program 

The  Ohio  State  University 
Princeton  University 

"Optically  Pumped  Nonequilibrium 

Plasmas" 

I.  Adamovich,  V.  Subramaniam,  W*  Rich 

Presented  at 

AFOSR  Workshop 
Understanding  and  Control  of 
Ionized  High-Speed  Flows” 

Princeton  University 
Feb.  26-27, 1998 

Work  funded  by  the  Director  of  Defense  Research  & 
Engineering  (DDR&E)  within  the  Air  Plasma  Ramparts 
MURI  Program  managed  by  AFOSR 


Nonequilibrium  Thermodynamics  Laboratories 


Program  Objectives: 


•  Create  large  volume  (0[m3])  free  air  plasma 

•  One  atmosphere 

•  Free  electron  densities  1013  cm3  or  greater 

•  Gas  temperature  less  than  2,000  K 

•  Energy  efficient  method 

Approaches: 

A.  Thermal  Plasmas  (High  Temperature  Arcs,  Reentry 
Plasmas) 

Problems: 

S)  High  Temperature 

*B.  NonThermal  Plasmas  (Glow  Discharges,  Optically 
Pumped  Plasmas) 

Problems: 

I)  Requires  thermodynamic  work  (electrical  work, 

laser  work) 

II)  Stability 

OSU/Princeton  Approach:  Laser  Pumped  Metastable 
State  Plasmas. 


Nonequilibrium  Thermodynamics  Laboratories 


The  Ohio  State  University 


Nonequilibrium  Thermodynamics  Laboratories 


The  Ohio  State  University 


Advantages: 


Stable,  large  volume 
High  pressure 
Low  gas  temperatures 

Problems: 

Electron  density 
Efficiency 
Non-air  species 


Nonequilibrium  Thermodynamics  Laboratories 


Three  Level  Laser 


{ 


Energy 

A 


Pumping  Waste  Heat 


Upper  level  energy:  Eu 
Lower  level  energy:  El 

Ground  level  energy:  Eo=0 


Quantum  Efficiency 
Ey  ~EL 

*1  quantum  g 

L 

~  40%  for  CO2/N2 


95%  forCO/N2 


Scoville,  1954  showed  equivalence 
with  reversible  engine  efficiency 


Poreur/AL  enc/ZGy  (fLecr/io/v  votr^J 


IR/UV  Radiation  Conversion 
in  Optically  Pumped  CO-Ar  Plasma 


Relative  population 


Nonequilibrium  Thermodynamics  Laboratories 


The  Ohio  State  University 


IR/Visible  Radiation  Conversion 
•  in  Optically  Pumped  CO-N2-Na  Plasma 


Relative  population 


Vibrational  Quantum  Number 


Nonequilibrium  Thermodynamics  Laboratories 


The  Ohio  State  University 


Nonequilibrium  Thermodynamics  Laboratories 


The  Ohio  State  University 


Nonequilibrium  Thermodynamics  Laboratories 


The  Ohio  State  University 


GRADIENT  DISTORTION  EXPERIMENT 

•  Transverse  temperature  gradient  is  created  by  absorption  of 
CO  laser  radiation  (1-2%  CO  is  added  to  the  flow  to  enable 
absorption  and  as  a  thermometric  element). 

At  P=100  Toit,  dT/dr  ~  100  K/cm 

Routinely  used  in  CO  and  NO  optical  pumping  kinetic 
experiments. 

•  Plane  shock  wave  created  by  a  capacitor  bank  discharge. 
Measurements  oft 

•  flow  density  distribution 

•  Gas  temperature  and  transverse  temperature  gradient  before 
the  shock  arrival 


Radial  temperature  profile  is  parabolic: 
T(r)  =  Tc  -  (Tw-Tc)(r/R)2 


Nonequilibrium  Thermodynamics  Laboratories 


Non-Equilibrium  Thermodynamics  Laboratories  &  Center  for  Advanced  Plasma  Engineering 


f 


Non-Equilibrium  Thermodynamics  Laboratories  &  Center  for  Advanced  Plasma  Engineering 
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